pos_proportion_mat = NaN(14,23);

for i = 1:23
    index = data(:,1) == i;
    n_subclass = sum(index);
    for j = 3:16
        n_pos = sum(index .* filtered_data(:,j));
        ratio = n_pos/n_subclass;
        pos_proportion_mat(j-2,i) = ratio;
    end
end

for i = 1:23
    base = i*ones(1,14);
    total = [total, base];
end


    